import numpy as np
import scipy.linalg as sla
A = np.array([
[2,1,1,0],
[4,3,3,1],
[8,7,9,5],
[6,7,9,8]
])
np.set_printoptions(precision=2, suppress=True)
print(A)
P, L, U = sla.lu(A)
print(P.T)
print(L)
print(U)
print(np.dot(P.T, A))
print(np.dot(L, U))